The Complete Jamming Landscape of Confined Hard Discs 
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An exact description of the complete jamming landscape is developed for a system of hard discs 
of diameter a, confined between two lines separated by a distance 1 + ^/3/4 < H/a < 2. By 
considering all possible local packing arrangements, the generalized ensemble partition function 
of jammed states is obtained using the transfer matrix method, which allows us to calculate the 
configurational entropy and the equation of state for the packings. Exploring the relationship 
between structural order and packing density, we find that the geometric frustration between local 
packing environments plays an important role in determining the density distribution of jammed 
states and that structural "randomness" is a non-monotonic function of packing density. Molecular 
dynamics simulations show that the properties of the equilibrium liquid are closely related to those 
of the landscape. 

PACS numbers: 64.70Pf, 61.43Fs, 64.60.My,63.50Lm,64.70qd 
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The properties of a wide variety of materials, includ- 
ing liquids, glasses, crystals and granular materials, de- 
pend on the way the particles pack. Bernal [l[ origi- 
nally used random packings of ball bearings to study the 
structure of liquids and coined the term random close 
packing (RCP) to describe the most dense random ar- 
rangement of spheres. Subsequently, the potential energy 
landscape (PEL) has become an important paradigm 
used to describe the role of particle packing in both the 
thermodynamics and dynamics of many of these sys- 
tems [H, Q] . I n this approach, each configuration of the 
liquid is represented as a point in the high-dimensional, 
TV-body, potential energy function of the system that can 
be uniquely mapped to the closest mechanically stable 
packing or inherent structure 0. For a system with a 
soft potential, an inherent structure represents a local po- 
tential energy minimum. For a hard particle system, an 
inherent structure is a collectively jammed packing [y, @| 
and a local density maximum, which gives rise to the 
corresponding jamming landscape (JL). The configura- 
tions that map to the same inherent structure can be 
grouped together into a local basin of attraction and the 
properties of the liquid can be described in terms of the 
number of inherent structures and the motion of the sys- 
tem through the resulting landscape of basins and saddle 
points. 

It has been suggested that RCP might be related to 
an ideal glass state. In the context of the JL, the config- 
urational entropy of a hard sphere system is defined as 
S c = lnf2(0j), where £l(<f>j) is the number of collectively 
jammed states with an occupied volume fraction <f>j, so 
the ideal glass transition corresponds to a density where 
the mctastablc liquid would become trapped in a single 
unique basin and S c — > 0. However, obtaining a detailed 
description of the JL remains a considerable challenge. 
Computer simulation has been used extensively to study 
packing, but different protocols often lead to different 
conclusions regarding the density distribution of inher- 
ent structures for both hard disc mixtures [H, 0, [l(| and 



hard spheres 11,0]. A recent study [l3[ of jammed 
packings also raised questions concerning the relation- 
ship between the structure of a packing and its density. 
In particular, it has been suggested that the RCP should 
be replaced by a maximally random jammed (MRJ) state 
that is more rigorously defined with respect to a set of 
order parameters. The only exact analytical results avail- 
able for the entire JL are for one-dimensional 1J| or 



quasi-one-dimensional [15[ hard particle systems where 
particles can only interact with a s ing le neighbor on each 
side. A recent mean field theory [1 61 ] . which finds hard 
sphere packings ranging from <pj = 0.536 — 0.635, repre- 
sents one of the few analytical results describing the JL 
in higher dimensions. 

The goal of the present work is to obtain an exact de- 
scription of the JL for a system so that we can explore 
the relationship between packing structure and density. 
To this end, we study a system of two-dimensional (2d) 
hard discs of diameter a, confined between two hard walls 
separated by a distance of 1 + \/3/4 < H/a < 2. In 2d, 
a particle is locally jammed if it has at least three rigid 
contacts that are not all in the same semicircle. How- 
ever, local jamming of all the particles is a necessary but 
not sufficient condition for collective jamming because 
the concerted motion of a number of particles can lead 
to a collapse of the structure Q . By confining the discs 
to a channel with H/a < 2, we prevent the particles from 
passing each other and eliminate the possibility of collec- 
tive rearrangements. The number of structures in which 
all the particles are locally jammed is then equal to the 
number of collectively jammed inherent structures. For 
H/a < 1 + iy/3/4, discs can only contact their nearest 
neighbors, which only allows two local particle arrange- 
ments that jam and gives rise to a binomial-type density 
distribution of inherent structures [HI . 

In the range of channel diameters studied here, the 
discs can contact both their first and second nearest 
neighbors, leading to a significant increase in the num- 
ber and structural diversity of possible locally jammed 



environments. The local packings are identified using a 
heuristic algorithm that generates new packings by mak- 
ing perturbative moves of individual particles, or groups 
of particles, on the known jammed structures for a system 
with H/a = 1 + \/3/4. Our rationale for the algorithm is 
based on the observation that the particle arrangements 
for 1 + \/3/4 < H/ a < 2 must collapse continuously back 
to one of the jammed structures at 1+^3/4 because there 
is no change in the number of possible contacts a disc can 
have over this range of channel diameters. 

All the local packing environments are then mapped 
onto a set of 32 tiles (Fig. fjji) that can be combined in 
a sequence, from left to right, such that the local jam- 
ming conditions for each particle can be achieved by just 
considering the neigboring tiles. This allows us to use 
the transfer matrix method to construct the exact parti- 
tion function for all the jammed states. Incompatibilities 
exist between some of the tiles in the sense that they can- 
not form a left-right pair that results in a valid, jammed 
packing of the particles (Fig. [TJa). This results in three 
groups of tiles; a set of high density tiles (tiles 1-5 in 
Fig. DJ<), so called because they appear in the high den- 
sity inherent structures; a set of low density tiles (tiles 
12-16) and a set of interface tiles (6-11). There is no 
direct compatibility between the low and high density 
tiles except that interface tiles are compatible with some 
members of both other sets, so it is possible to mix the 
groups within a single packing. The separation in the 
packing compatibility arises because the two important 
length scales in the system, a and H, are incongruent. 
Tiles 12 — 16 all contain a disc-disc contact that spans 
the width of the channel at an angle that cannot jam the 
particles in the other set of tiles. 

We begin to construct the partition function by defin- 
ing the 32 x 32 transfer matrix, M, with matrix elements 
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where = 1 if tiles i and j are compatible but is zero 
otherwise. n.y and are the number of particles and 
length associated with the addition of tile j to the right 
of tile i. [i is the chemical potential and Pls is the exter- 
nal longitudinal landscape pressure applied to the ends 
of the channel, rather than the pressure related to the 
internal vibration of the particles which must be zero for 
the jammed states. f3 = 1/kT, where T is the temper- 
ature and k Boltzmann's constant. The definition of a 
temperature for a system of jammed states is of consider- 
able general interest However, here we note that the 
configurational integral of a hard particle system is inde- 
pendent of T so we are free to simply treat temperature 
as a parameter conjugate to the entropy. Furthermore, 
an analysis of the transfer matrix shows that the num- 
ber of states with a given N and V appear as coefficients 
of the polynomial matrix elements of M Nt that can be 
extracted by taking the appropriate derivative with re- 




6 J^__ 8 9 10 11 



16 



12 13 14 15 




FIG. 1: (a) The set of tiles that represent the local jam packed 
configurations. The vertices represent disc centres and the 
solid lines join the centres of two contacting discs. The small 
open circles identify those discs that contact the wall. An 
additional 16 tiles are generated by reflecting each tile in the 
plane of symmetry that runs through the axial center line 
of the channel, e.g. tile 1' is generated from tile 1. (b) A 
configuration constructed from tiles 2' — 1 — 12. A tile is 
added to the sequence by sharing the first two discs with the 
previous tile, with the exception of tile 14, which only shares 
one disc, and tiles 6-8, which may either simply contact the 
previous tile, or share discs, depending on the nature of the 
neighbouring tile. In our example, tile 1 is compatible with 
2' (CV.i = 1) but tile 12 is incompatible with 1 because it 
causes particle overlap between the unfilled discs (Ci,i2 = 0). 
Other incompatibilities result in discs remaining unjammed. 
(c) Most dense packing. 



spect to T without necessarily providing a clear thermo- 
dynamic definition of temperature. 

Taking M Nt gives all the possible jammed packing 
arrangements that can be formed with Nt tiles. Both 
the number of particles in a packing, N, and its volume, 
V, fluctuate between different packings, so we use the 
generalized ensemble partition function 18| . which can 
be expressed 



r(p LS ,fx,T) 
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where A^ is the i th eigenvalue of M. The equilibrium 
condition is obtained by finding P£ s and fj,* such that 
T(P£ S , fx*,T) = 1. In the thermodynamic limit, the 
equilibrium properties of the ensemble averages are over- 
whelmingly determined by packings of a given N and V, 
which arc obtained from. 



N = c^nr^,^) | ( p, s ,m)->(p^) (3) 
V = dp LS lnT(P LS ,v,T)\ {PLSt ^ (P . ,„.) , (4) 
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FIG. 2: Configurational entropy, S c vs cj>j. Insert: Enlarged 
high density region. 

and the configurational entropy, 

S c = kTd T hiT(P L s,fi,T)\ { p LS ^ (P * s ^, ) , (5) 

is the logarithm of the number of such packings. We 
also calculate the fraction of particles, /y, belonging 
to interfaces between tiles i and j, as follows: if an 
additional factor 7 is associated with matrix element 
Mij = expi-pPLsHLij + npjNij) then, 

fij = (V 'nN)d 7 lnT(P LS , /j,,T,'y)\(p LStll) _>(p* s!tl *y (6) 

All results reported here are given for a channel width 
of H = 1.95cr. Fig. [2] shows S c /Nk, which is indepen- 
dent of T, as a function of the <fij over the entire range 
of packing densities, except at the very extremes where 
it was not numerically possible to find solutions where 
Pls and A* becoming extremely large. Both the most 
dense and least dense states are non-degenerate, except 
with respect to identical structures obtained by a reflec- 
tion along the axis of the channel. The most dense state, 
with (f>j f=a 0.80743, consists of a repeated tile sequence of 
— 1 — 2 — 1' — 2'— (Fig. [Tt), while the least dense pack- 
ing is made from alternating —15 — 15'— tiles, giving 
(f)j ks 0.613. The maximum in S c occurs at <f>j = 0.712, 
which also corresponds to the density sampled by the sys- 
tem at Pls = (see equation of state (EOS) in Fig. [5k.). 
This is consistent with earlier studies that found config- 
urations of the ideal gas were related to the maximum in 
the configurational entropy. The lower density packings 
are sampled at negative pressures, which are non phys- 
ical in an equilibrium hard particle fluid, while the sys- 
tem samples deeper basins with increasing Pls- Above 
4>j = 0.806, we see a smooth, but rapid change in the 
EOS as S c dramatically decreases. 

The structure of the packings is examined by measur- 
ing the fraction of particles associated with the high den- 
sity tiles, fhdi low density tiles, fid-, and interface tiles, 




FIG. 3: (a) <j>j as a function of Pls- (b) <j>j as a function 
of the equilibrium fluid <j> obtained from simulation studies 
with N — 996 (circles) and N = 7998 (squares) particles, 
(c) The fractions of particles in tile interfaces between high 
density tiles, fhd = Yl^j—i fa ( sou d line), low density tiles, 
fid = j-12 fa (dashed line) and interface tiles, fmt = 
X^i=6 6 *=r fa (dot-dashed line), (d) high density enlargement 
of (c)."~ 



fint, as a function of <j>j (Figs. [5fc and [3p!). The sys- 
tem can increase <f>j in two ways: by replacing the more 
open tiles of a set with more dense tiles from the same 
set, i.e. replacing tile 15 with tile 9; or by replacing low 
density tiles with high density tiles. Fig. shows that 
both these processes are occurring over the entire density 
range, but we also notice that fi„t is always decreasing 
above <pj = 0.64. This implies that the packings are 
becoming more structurally heterogeneous and contain 
clusters of high and low density tiles as their ability to 
mix is reduced. By <pj = 0.806, most of the low density 
tiles remaining are of type 9, so that the only way to 
increase the density further is to eliminate the low den- 
sity tiles altogether, which results in a rapid structural 
crossover and the decrease in S c . 

Structural randomness would imply that it is equally 
probable to find any given packing arrangement, on 
all possible length scales. Consequently, the function, 
* = (l/d ma x)'E dt i t jPij(d)'ln-Pij(d), where Pij(d) is the 
probability of finding tile j a distance d from tile i and 
the sums are over all possible tile pairs and separations, 
should be a maximum for the most random states. d max 
is the maximum separation studied. ^ is a measure of 
the average randomness of the packings at a given <f>j 
and is likely to provide a good description of the liquid 
structure as a function of density. In contrast, the MRJ 
state focuses on the randomness of a particular structure. 
Fig. d] shows the unexpected result that randomness is 
not a monotonic function of <j)j. The global maximum 
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FIG. 4: * as a function of fa obtained at T = 0.15. The 
diamond represent * for the most dense state, obtained di- 
rectly from its structure. Insert: An enlargement of the high 
density region. 



in 'F occurs at <f>j ps 0.66 and it initially decreases with 
increasing (f>j as a result of the clustering of low and high 
density tiles described above. The elimination of the low 
density tiles then allows greater mixing between the high 
density tiles, causing 'F to increase again to its second 
maximum at <f>j ~ 0.8071. 

We use event driven, molecular dynamics (MD) sim- 
ulations [psl ] to understand how the thermodynamics 
of the equilibrium fluid are related to the landscape. 
The system is decompressed in small increments from 
<t> = 0.8072, using the most dense packing as the starting 
configuration, to <f> = 0.1. At each 4> studied, 100000 N 
collisions are used to reach equilibrium and the data are 
collected over the next 25000A?" collisions, where N is the 
number of discs in the simulation. 200 independent con- 
figurations from each </> are then compressed until they 
are jammed using a modified version of the Lubachevsky 
and Stillinger [20| (LS) algorithm that ensures H/a re- 
mains constant as the diameter of the discs is changed (L 
fixed). A compression rate da/dt = 0.04 is used. This 
process is not an "infinitely" fast quench and the sys- 
tem is able to relax via MD while it is compressed, but 
faster quench rates result in poorly jammed states with 
loose discs. Fig. |3Jd shows that the average 4>j, obtained 
by quenching the equilibrium liquid from <j) for systems 
with N = 996 and A" = 7998 discs, displays a dramatic 
change in its 4> dependence above (f> ss 0.66. This change 
is also apparent in the EOS (Fig. [3^) for the jammed sys- 
tem obtained directly from the partition function. While 
the EOS always varies smoothly, signifying that there is 
just one thermodynamic fluid phase at all densities, the 
rapid change in slope of <j>j coincides with the structural 
crossover observed in Fig. [3)1. 

More generally, our work highlights the importance 



of understanding the role of local packing environments 
and geometric frustration in determining the properties 
of hard particle packings, and its relationship to pack- 
ing randomness. Anikeenko et al. [13, 21 1 also observe 
competition between the "quasiperfect" tetrahedra of the 
disordered packings and the crystalline arrangements of 
hard spheres. This suggests frustration may be a generic 
feature of packings and that a tiling approach may be 

We also 



22] 



useful in enumerating inherent structures 
find that S c is exponential in the number of particles over 
the entire density range of (f>j, which is consistent with 
Donev et al. [lfj, who used computer simulation show 
that this is the case for binary mixtures of hard discs in 
the bulk. Finally, the complete knowledge of the JL for 
a system, along with an understanding of local packing, 
will allow us to directly test the relationships between the 
landscape, kinetically facilitated dynamics geomet- 
ric frustration [24|, and the glassy dynamics of liquids. 
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